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Abstract 

An extended version of the Nambu-Jona-Lasinio (NJL) model is applied to describe 
both nuclear matter and surface properties of finite nuclei. Several parameter sets are 
discussed and a comparison of the saturation properties and equation of state (EOS) 
with the NL3 parametrization of the non-linear Walecka model is made. The properties 
of asymmetric matter are discussed. 



1 Introduction 



Strong interaction dynamics of mesons and baryons is believed to be described by QCD. This 
theory exhibits a non perturbative behaviour at low energies, a circumstance which renders 
the analytic study of the theory rather difficult. The NJL model is a popular substitute 
which has in common with QCD important symmetries of the quark-flavour dynamics. This 
model has been very successful in the description of the meson sector. The question arises 
whether the model also allows for inhomogeneous solutions at non-zero density which are 
of interest for the description of hadronic matter. 

The NJL model ^ was originally developed for the purpose of understanding hadron 
physics. In this model, hadron masses are generated by spontaneous symmetry breaking 
of the vacuum. In the modern form of the NJL model, we start with essentially massless 
quarks interacting via zero range interactions, but with a cutoff in momentum space [21 1313 • 
NJL is thus only an effective theory, in which form factors and other finite range effects 
have been ignored, and it does not take into account quark confinement. However, it may 
work quite well in the region of interest in nuclear physics, i.e. for excitations less than the 
scalar meson mass. 



The paper has been organized as follows: It is shown that with the generalization of the 
NJL model for nucleons, referred to as the Extended NJL (ENJL) model, it is possible to 
get a reasonable nuclear equation of state and behaviour of the effective mass. Then we 
list some numerical results, concerning both nuclear matter and quark matter. It is found 
that the performance of the ENJL model is almost perfect as far as bulk static properties, 
such as saturation and binding energy of nuclear matter, or the effective mass of nucleons, 
are concerned. However, the model fails when we wish to deal with dynamical properties, 
since it leads to scalar excitations with the unacceptable mass of twice the nucleon mass, 
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or to a poor description of the nuclear surface. In order to overcome these drawbacks, the 
ENJL model is replaced by a chiral invariant model of nucleons interacting through the 
exchange of a , tt and oj fields, which we call extended chiral sigma (ECS) model. For the 
description of bulk static properties, the ECS model is almost equivalent to the ENJL model. 
In these models, the nucleons are regarded as composite particles. The composite nature 
of the particles reflects itself in the fact that the coupling depends on the local density. 
The surface properties of finite nuclei are calculated within the ECS model. Some remarks 
concerning the connection of the ENJL and ECS models to other relativistic chiral models 
El are presented. Then, some properties of asymmetric nuclear matter are discussed 
and some brief conclusions are drawn. Next, we discuss a two dimensional version of the 
NJL model for which some simple analytic results can be obtained. This is followed by a 
comparison of our model with several other models, including the quark meson coupling 
model and versions of relativistic mean field models. Preliminary results of the present work 
have been presented in [0]. 



2 NJL model 



2.1 Original NJL Model 



The NJL model |3j is defined by the Lagrangian density 

C = i^ii^d^,)^ + G[{4>^f + [i^ii^T^f]. (1) 

A regularizing momentum cut-off A is part of the model. The Lagrangian is equivalent to 
the Hamiltonian 

N N 

Hnjl = Y.Pk-ak + G 5] - - ihin ■ n) ■ (2) 

k=l k,l=l 

The vacuum is described by a Slater Determinant |<I>o) constructed from plane waves which 
are negative energy eigenfunctions of the single particle Hamiltonian h = p ■ a + (5m. The 
"constituent mass" m is a variational parameter. 



If moreover positive energy eigenfunctions with momentum p satisfying \p\ < Pf are oc- 
cupied, so that Pf is the Fermi momentum, the energy expectation value E = ("S^ol^Af jlI^o) 
is given by 



For quark matter, the degeneracy is = 2N(.Nf and A is such that m = 313 MeV is the 
constituent quark mass in the vacuum. The gap equation, which determines m, reads, 

E / ^ , =1- (4) 
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2.2 Nuclear Matter vs Quark Matter Cutoff 



It is well known that the NJL model does not contain a mechanism for quark confinement. 
Thus the clustering of quarks into nucleons is not accounted for in this model. In this paper, 
quark clustering is put in by hand. The way it is done here is that we treat the nucleons 
as the elementary particles, but take quark structure into account in the interactions. We 
will find it convenient to define dimensionless integrals In by: 

-dp = In{x) . (5) 



According to the NJL model, the volume integral of the A^A^ interaction is given by 

^' = = -S-f^^e'Mo^^l^)' . (6) 

Here N,, is the number of colors, and Mg is the free nucleon mass. For quark matter, N^. = 3. 
The scalar (cr) meson mass is: 

^' = '^^ = ^- 

Let us now try to get the same A^A^ volume integral with a single color. This can be 
accomplished by using a lower value of the cutoff, denoted here by A'. Further, from here 
on, we will often express all masses in units of the free nucleon mass Mq. Also, we will set 
h and c equal to 1. We find that: 

,,,AO = M^. (8) 

An alternative way of taking quark clustering into account is by fixing A ' so that the free 
nucleon mass in vacuum, Mq is preserved. Actually, such a prescription has been adopted 
in some of the calculations which are reported in Section 4 of this paper. However, the value 
of A ' so obtained is in a very reasonable qualitative agreement with the one determined by 
eq. ©. 



2.3 Extended NJL Model 



The NJL model can be extended |71 |H1 El El to yield reasonable saturation properties 
of nuclear matter, the field ip being then the nucleon field. In Ref. 0, the nucleon is 
constructed as a 3-quark bound state, this quark substructure giving rise to a mechanism 
for saturation of nuclear matter which plays the same role as the 8-fermion term in our 
ENJL model. Such substructure may therefore be regarded to provide a microscopic basis 
for our phenomenological approach. An effective density dependent coupling constant is 
obtained if the following extended NJL (ENJL) Lagrangian density, which actually pushes 
chiral symmetry restoration to higher densities, is considered, 

- Gsv[{i>il^f + {i>il^T^l:)\i>i^i^f. (9) 
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Coupling constants, etc. 



Quantity 

G,(fm2) 
G,(fm2) 
G..(fm8) 
Gp(fm2) 



EOSI 

3.880 
3.952 
-4.901 
2.794 
0.913 
0.427 
-1705.4 



h 

p/pQ (M = 0) 



Quantity 

Mo (MeV) 
A (MeV) 



EOSI 

939 
418.9 
0.0560 
0.0502 

3.0 

111 

4.3 



Esq = Ebn (MeV) 



Vs 



p/po {Esq = Ebn) 



Table 1: The coupling constants and properties for our preferred equation of state 

For nuclear matter "0 is the nucleon field, the degeneracy \s u = 2Nf and A is such that 
M = 939 MeV is the nucleon mass in the vacuum, determined variationally, as explained 
below. A different extension which also leads to nuclear saturation has been made jlll I12j . 
This involves the fourth power of the scalar density. Another alternative way of extending 
the NJL model was studied by Rezaeian and Pirner |13j . Here the nucleon is considered to 
be a composite particle, a bound state of a quark and a diquark. The term in Gy is supposed 
to simulate a chiral invariant short range repulsion between nucleons connected, possibly, 
with the chromomagnetic interaction between quarks. The term in Gsv accounts for the 
density dependence of the scalar coupling between nucleons which is a manifestation of 
the composite nature of these particles and is equivalent to allowing the coupling constant 
G in the original NJL model to depend on p in such a way that G becomes replaced 
hy Gs{p) = Gs{l - %fp'^). It is found that the magnitude of the a field decreases as 
the hadronic density increases, in accordance with the corresponding tendency for chiral 
symmetry restoration. This mechanism has a similar effect, for saturation, as the analogous 
one proposed by Guichon according to which confinement is affected by an increase 
of the density in such a way as to become gradually less effective, which amounts to an 
increasing repulsion. For nuclear matter, the NJL model leads to binding but the binding 
energy per particle does not have a minimum except at a rather high density where the 
nucleon mass is small or vanishing. Thus the pure NJL model does not give a proper 
equation of state. The introduction of the Gsv coupling term is required to correct this. 
On the other hand, for quark matter we set G^ = Gsv = 0. For quark matter the energy 
per particle W predicted by NJL does not lead to binding and does not have a minimum 
except for unreasonably small values of the cutoff momentum A. 

The thermodynamical potential per volume corresponding to Q is 



Here M is the p dependent nucleon mass. To avoid a cumbersome notation, the p depen- 
dency is not explicitly indicated but is implicitly assumed, without danger of confusion. 
The free nucleon mass Mq is the value of M for ^ = 0. The condition duj/dpp = implies 




Ep^=p- 2Gv{^H) - 2Gsv{i^H){m\ 
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EOSI 


NL3 


E/A- Mn(MeV) 


-16.12 


-16.3 


Pn (fm~^) 


0.148 


0.148 


M/Mo 


0.75 


0.60 


K (MeV) 


295 


272 


Wsurf (MeV) 


19.43 


19.38 


R fm 


5.33 


5.328 


t fm 


2.64 


2.65 



Table 2: Nuclear matter saturation properties and surface properties 



with Epp = y M'^ These conditions fix the values oi pp, M for given fi. We observe 

that the gap equation, which determines M, is the same as eq. (0} except for the replacement 
of G with Gs{l- ^fp^)- The choice with Gs as a phenomenological parameter fixed so that 
nuclear matter saturation properties are reproduced, we will refer to as model I. The value 
of Gs is about 20% smaller than that obtained from the NJL model. On the other hand, it 
is also tempting to choose Gg for nuclear matter nine times bigger than the quark matter 
value, G. The philosophy behind the last assumption is that the NN interaction is due to 
the instanton interaction between quarks, predicted by QCD in the weak coupling regime. 
This choice of Gs defines what we will denote in the sequel as model II. 

The properties of the extended NJL model are now easily computed. 



3 Properties of nuclear matter 



In Table numerical results pertaining to our preferred model EOSI are presented. As 
input, we use the saturation properties of nuclear matter fitted empirically in Ref.|15j. 
The three coupling constants Gs-, Gy and Gsv are adjusted to fit an effective mass M{po) = 
0.75Mo and the saturation energy and density. Model EOSI seems to be quite realistic, since 
it fits well both saturation properties and the nuclear surface properties. We also compare 
our results with those of T![^ in which a relativistic mean field theory (RMF) is used, and 
which describes well both stable and unstable nuclei. The RMF parametrization introduced 
in is known as NL3. Model I gives results quite close to NL3 especially for the surface 
properties (see Table I^J. The relation between the coupling constants Gs, G^, Gsv and the 
parameters c^, Csv, which will be used later, is fixed as follows: 

^ 2/i(A0M2 ^^^^ 
G - \ ^ C (12) 

6 4 

f-l ^ ^sv _ _p ^ ^SV /-.o\ 

'^^'^ - 2/i(A')2/2(A')Af8 '^^/i(A')/2(A')M6- ^'^^ 

Here Mq = 4.76/m~^. The crucial role played by clustering must be stressed. It should 
be pointed out that without clustering there is no real binding and that, moreover, the 
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P/P„ 

Figure 1: Energy vs density in ENJL Model for symmetric nuclear matter (EOSI and 
EOSII) and for quark matter (EOSI-Q) 

incompressibility becomes unacceptable. In Fig. ^we plot the binding energy for nuclear 
matter and quark matter (EOS-Q). For the latter case, we use a quark NJL model which 
differs from the nucleon NJL model (ENJL) used in the rest of this paper. At some density, 
the nucleon matter curve intercepts the quark matter curve, so it is clear that, at high 
density, quark matter prevails and may be found, for instance, in the core of neutron stars. 
Our experience with this problem tells us that without clustering and without the eight 
fermion interaction (that is, the density dependent coupling) saturation comes closely after 
chiral symmetry restoration, almost coinciding with it. Clustering and the eight fermion 
interaction places chiral symmetry restoration well after saturation. 

The philosophy underlying model EOS II is to use a coupling Gg 9 times bigger than G in 
the quark NJL model and to fit one less quantity. Thus the effective mass at saturation 
was left free. We have taken E/A = -15.8 MeV and po = 0.148 fm"^, i.e. kp = 0.130 fm~^ 
according to jl5l I16j . Using this model, one gets quite different results for the coupling 
constants from EOS I, namely Gs = 1.746 fm^, Gy = 3.387 fm^, Gsv = -1-839 fm^, A = 553 
MeV . The effective mass at the saturation density turns out to be 0.89, much closer to the 
free nucleon mass than for model I. Also, this model leads to surface energy and thickness 
only about half of the empirical values. Similar results are obtained by Koch et al, ref. jjj. 
However, the cutoff A obtained using EOS II is not far from the one determined by eq. (jSJ. 



4 Surface properties 

Next we will calculate nuclear surface properties which are mainly determined by the scalar 
meson, the meson responsible for the nuclear attraction. We have seen that the ENJL model 
accounts adequately for bulk static properties of hadronic matter, such as saturation and 
binding energy of nuclear matter, or the effective mass of nucleons. However, the model fails 
when we wish to deal with dynamical properties, since, as it has been shown (HI El) it leads 
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Figure 2: Effective Mass vs density in ENJL Model for models I and II and results of Koch 
et al. [7] 



to scalar excitations with the unacceptable mass of twice the nucleon mass, and to a poor 
description of the nuclear surface. We keep the ENJL Lagrangian, in spite of the referred 
drawbacks, because it is extremely simple and transparent, clarifies the underlying physics 
and incorporates in the simplest possible manner the basic philosophy of our approach. 
In the extended chiral sigma model, introduced in ^ ^], the nucleons and mesons are 
regarded as composite particles. The version used in J2] is referred to as ESM. In the 
present paper and in |^, a slightly different version is used which we refer to as ECS. Both 
models are equivalent to corresponding ENJL models, as far as ground state properties are 
concerned, provided exchange terms are neglected. However, there are significant differences 
between the models. In ESM, the coupling to the sigma meson depends on the scalar- 
pseudoscalar field, while in the ECS model, that coupling depends on the baryon density. 
Moreover, the nuclear calculations are more easily implemented in the ECS model, which 
is free of instability that may occur in ESM. The advantage of the ECS model over ENJL 
stems from the fact that it generates adequate collective properties, this being the reason 
that ECS can be used to study the nuclear surface. 
We write formally the ECS model Lagrangian density as 



_ _ _ _ _ 1/2 ~ 



with 



2mj zmf, 2mj 

For comparison, the ESM Lagrangian density is given by Eq.(ll) of [T21, while in the linear 
sigma model, ref. ^7], we have: 



1 

2 

^2 + ^ . 7? - F^f (16) 
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Here the nucleon mass is generated by a "Mexican Hat" potential. If we write a = aQ + F.,^ 
then the last term in the linear sigma model contributes — |m^crQ(l + ^-)^ to the La- 
grangian. It was shown in jl8j that the Dirac-sea in the linear sigma model with valence and 
Dirac-sea quarks but no "Mexican hat" provides the above effective mesonic self-interaction. 
The quantity -F^ is equal to o" = M/g^, when the energy is minimized w.r.t. a. For zero 
particle density, we obtain, according to the ECS model, Eq. ((Til) . 

9s ^9s Q 



2ml 8M2 ^ 2/i(4)M2 

which leads to 

F^ECS) 



3Mn 



A iV2 



(17) 



2tt 

This is to be compared to the pion decay constant according to the ENJL model, which, 
for vacuum properties such as the present one, coincides with the nucleon NJL model, 



FANJL) = ^ 
z vr 



1/2 



(18) 



We observe that Eq. (fTH)) may be obtained from Eq. (3.26) (a) of Ref. jSj upon replacing the 
number of colours Nc by 1 and m* by Mq. The quantities Ii and I2 are defined in Section 
2.2. The calculated pion decay constants for the ECS model (106.1 MeV) is in fairly good 
agreement with the empirical value of 93 MeV, but the pion decay constants for the ENJL 
model (33.49 MeV) is not. This was expected since the performance of the ENJL model is 
poor for dynamical properties. 

In order to describe surface properties we must include gradient terms in the meson field 
[T^ . Somewhat analogously to our approach, in Ref. supplementary gradient terms 
were introduced into some kind of extended NJL model with nucleons. This is in consonance 
with the work 001 where an effective field theory of many-nucleon systems is obtained from 
the NJL type Lagrangian. In the ECS Lagrangian ()14() . we will keep only the gradient 
terms corresponding to the scalar meson, the philosophy behind this approximation being 
that the relevant length scale is determined by the sigma mass. 

The corresponding thermodynamical potential per volume in the Thomas-Fermi approx- 
imation is then 

ioiii) = (^(7.p»+5,a(^V)(l-^(V'^V'>'j (19) 



Minimization with respect to a, Vq yields eq. (|1U|) for infinite matter, provided we choose 
Gs = g1/{2m1). Variation with respect to M, pp-, c and Vq gives: 

M = 5.a(l-^(V'W)'^', 

Ep, = ^i- 5. Vo - gscj^m (V'V) (1 - ^{^^^?) > 
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E 




Figure 3: Density profile for nuclei with A = 100 



= m„Vb 



-gs{i'i>)[i 



1/2 



(20) 



Next, we look for a droplet solution within model I and compare its surface properties 
with the results obtained with the NL3 parametrization JS]. In the small surface thickness 
approximation (V^cr 



the free energy F can be rewritten as j21j : 



(21) 



where C is a constant and A is the number of particles. For droplets with radius R and 
volume y, 

F{R) = Wsurf A"^/^ -CV + fxA. (22) 
The surface energy per unit area of these droplets in the small thickness approximation is 



surf 



(23) 



The surface thickness t is defined as the width of the region where the density drops 
from 0.9pBo to O.lpsoi where pso is the baryonic density at r = 0. In Table [21 the val- 
ues of Wsurf ^ t, for a nucleus with A = 100, predicted by model I is displayed and 
compared with the corresponding results obtained by the NL3 parametrization. In the 
present approach the mass of the a meson appears as an extra parameter. We have chosen 
nis = 2mq = 626 MeV, as in the NJL model. In Fig. EJwe plot the density profile of a 
nucleus with 100 particles. 



5 Isospin asymmetric nuclear matter 



The recent advance in unstable nuclear-beam experiments has been providing us with new 
knowledge on unstable nuclei away from the stability line. We hope to get information on 
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the properties of dense matter under extreme conditions as the neutron-rich enviroment in 
neutron stars. 

In order to describe isospin asymmetric nuclear matter within the present formahsm it 
is important to include in the Lagrangian density the isovector-vector term 



■^p — ~Gp 



Here, Gp is chosen in such a way that the experimental value of the symmetry energy 
flsj/m = 35 MeV is reproduced. 

In Fig. |l]we compare EOSI with two different parametrizations of the non-linear Walecka 
model, NL3 JHI and TMl j22j . and the quark meson coupling model (QMC) 14 for sym- 
metric matter. In QMC the structure of nucleons is taken into account having as underlying 
model the MIT Bag model. In Fig. [21 we make a similar comparison with NL3 and TMl 
for neutron matter. The RMF parametrizations were fitted to both stable and unstable 
nuclei. However TMl contains a self-coupling term in the vector meson which weakens the 
contribution of the vector meson at high densities. We conclude that EOSI is softer than 
NL3. Indeed, EOSI is also softer than TMl and QMC, but becomes slightly stiffer at higher 
densities, in connection with chiral symmetry restoration. 

The nuclear EOS, as function of density and asymmetry, is reasonably well approximated 

by: 

^ = W=m[{-2{l-5)-^ + i^f] (24) 
where 5 = {pn — Pp)/{pn + Pp)- Here, pi, i = n,p is the density of particles of type i. 



6 Original NJL Model in Two Dimensions 

For the two dimensional NJL model, we get simple analytic results. The average kinetic 
energy T is given by: 
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p /p 



Figure 5: Binding energy versus density according to EOSI, NL3 and TMl, for neutron 
matter 

T = Pf — ^ p'j + . . ., where pf = 2tt pn- Thus we get, {cgv = for pure NJL): 



m 



W 



1 - 9Pf 
9-1 



Pf + Op} + 



(25) 
(26) 



where 5 = and m is the quark masss. For the two dimensional case, we can obtain exact 
results for these two quantities, not just for low densities: 



m 



W 



l-2gpf + 

9-1 
^P/- 



(27) 
(28) 



It is remarkable that higher order terms in W, cancel exactly. The vanishing term of order 
/9j in Eq. I\28\i has three components which cancel exactly: 

1. Relativistic correction to energy momentum relation. 

2. Increased kinetic energy due to reduced effective nucleon mass. 

3. Increased attraction due to reduction in effective mass of scalar meson (similar to that 
of the nucleon). 

For the more realistic three dimensional case, the physics is very similar, and again there 
is a near cancellation between these three terms. It cannot, of course, be perfect, since 
in three dimensions, these terms are proportional to and p'^, respectively. This 

explains why the Cgv term is essential for achieving saturation. 



7 Connection with other Models 



7.1 Guichon quark meson coupling model 



One promising possibility is the quark-meson coupling model proposed by Guichon jl4j . 
In this model, the nucleon mass is generated dynamically using a MIT bag model. Let us 
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define the dimensional coupfing constant: 



'~ mlM~ M ' 

where is the volume integral of the A^A^ interaction, eq. Then in the Guichon model, 
we automatically obtain a saturating term which can be put into our form: 

Wp2 = Csv BI M (? (30) 

where Cgv = 0.340 r\,ag , and P = The coefficient 0.340 is a dimensionless constant 
which can be derived from the MIT bag model . In that work, nuclear matter saturation 
and some properties of finite nuclei are fitted with a bag radius of r^ag = 0.8/m. This means 
that: Csv = 0.340 • f^f • 0.8 = 0.432, essentially the same as 0.427 obtained using the ENJL 
model. Their value for the volume integral of the scalar coupling is = — 1719 MeVfm^, 
very close to the ENJL result, Vs = -1705 MeV fm^. 



7.2 Relativistic Mean Field Models for Equation of State 

In the non-linear and derivative coupling mean field models, we also get a term quadratic 
in the density. In the strong coupling limit, neglecting kinetic energy, we can write for the 
energy per particle: 

1 fM-l)2 

The effective mass is obtained by setting ^ = 0. 

In the low density limit, the effective mass is given by: 

M = 1 - Bsp + ^aB^^p^ + ... (32) 

and 

W = ^{Bv- Bs)p+ \aBlp^ + ... (33) 

Making a comparison with Eq. ()3U() . we obtain Cgv = ■^a. In the original derivative 
coupling model we get a = 2, and Csv = 1, which seems too large. However, in the 
hybrid model developed by Glendenning et al [2^, we have a = 1, which corresponds to 
Csv = 0.5. This is close to our phenomenological result. On the other hand, with the other 
EOS parametrizations, Cgv has different values. 



7.3 Connection of ECS model to other relativistic chiral models 

In the present form the Lagrangian density (|14() contains terms coupling the meson fields 
to the fermion fields in a non-linear way. Using the equations of motion and keeping only 
quartic terms in the meson fields we get. 



+ \^-^a,{a^ + ^-^)V^^V^. (34) 
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This Lagrangian density, except for the meson kinetic terms, is similar to the one intro- 
duced by Boguta which, by including a scalar-vector coupling, was able to reproduce 
nuclear matter saturation properties. However, in an effective mesonic self-interaction 
corresponding to the "Mexican hat" is present. 

In a generalization of the model proposed in [H], including a "bare" vector mass, 
a quartic term in vector field and allowing for different vector-scalar and vector-nucleon 
couplings, was studied. It was argued that a linear realization of chiral symmetry was too 
restrictive and, although nuclear matter properties may be reproduced, there are problems 
with the properties of finite nuclei. In our formulation, however, and with our choice of 
parameters we have shown that not only nuclear matter saturation properties are reasonably 
explained, but also surface properties. 

8 Conclusions 

We have studied an extension of the NJL model which yields reasonable saturation of nuclear 
matter. An effective density dependent coupling constant is obtained which pushes chiral 
symmetry restoration to higher densities. Two variants of the model defined by appropriate 
sets of the parameters, Gs, G^, Gsv, and A, have been studied numerically. For Model I, the 
values of the defining parameters, the vacuum properties, the properties at saturation are 
given in Tables ^ and l2j An extension of the chiral sigma model has been developed which 
is based on the extended NJL model and fits not only the properties at saturation but also 
the empirical nuclear surface energy and thickness. The term in Gsv, responsible for the 
density dependence of the effective coupling constant, plays an important role in pushing 
to higher densities the restoration of chiral symmetry and in lowering the incompressibility. 

The relation of the present model to the relativistic chiral model proposed by Boguta [H] 
was presented. Finally we have discussed asymmetric matter within the present model. 

We have shown the connection of the ENJL model with other relativistic models including 
QMC, ECS, non-linear and derivative coupling mean field models. In order to further test 
the model it is important to calculate other properties of finite nuclei and its performance at 
finite temperature. The critical temperature for chiral symmetry restoration at zero density 
is easily calculated to be Tc = 196 MeV, an encouraging result. Application of the model 
to neutron-star matter requires the inclusion of strangeness and beta-equilibrium. 
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